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We develop a theory for two-dimensional diluted magnetic semiconductor systems (e.g. 
Gai-jMiijAs layers) where the itinerant carriers mediating the ferromagnetic interaction between 
the impurity local moments, as well as the local moments themselves, are confined in a two- 
dimensional layer. The theory includes exact spatial disorder effects associated with the random 
local moment positions within a disordered RKKY lattice field theory description. We predict the 
ferromagnetic transition temperature (T c ) as well as the nature of the spontaneous magnetization. 
The theory includes disorder and finite carrier mean free path effects as well as the important cor- 
rection arising from the finite temperature RKKY interaction, finding a strong density dependence 
of T c in contrast to the simple virtual crystal approximation. 
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Many projected applications of diluted magnetic semi- 
conductors (DMS), i.e. systems which combine the ad- 
vantages of a ferromagnetic material with those of a 
semiconductor with the additional flexibility of carrier- 
mediated fcrromagnetism enabling the tuning of the mag- 
netic properties by applying external gate voltages or op- 
tical pulses to control the carrier density, would involve 
the use of two-dimensional (2D) DMS structures such as 
quantum wells, multilayers, superlattices, or heterostruc- 
tures. Such 2D DMS structures are also of intrinsic fun- 
damental interest since magnetic properties in two di- 
mensions are expected to be substantially different from 
the three dimensional (3D) systems [1] that have mostly 
been theoretically studied in the DMS literature. The 
2D DMS systems introduce the possibility of gating, and 
thereby controlling both electrical and magnetic proper- 
ties by tuning the carrier density. In fact, such a car- 
rier density modulation of DMS properties has already 
been demonstrated [2,3] in gated DMS field effect het- 
erostructures. For various future spintronic applications 
the development of such 2D DMS structures is obviously 
of great importance. 

In this Letter, we provide the basic theoretical pic- 
ture underlying 2D DMS ferromagnetism focusing on the 
well-studied Gai-^Mn^As, with x ~ 0.03—0.08, a system 
where the ferromagnetism is well-established to be arising 
from the alignment (for T < T c ~ 100K - 200K) of Mn 
local moments through the indirect exchange interaction 
carried by itinerant holes in the GaAs valence (or impu- 
rity) band (that are also contributed by the Mn atoms 
which serve the dual purpose of being the impurity local 
moments as well as the dopants) . Our theory is quite gen- 
eral and should apply to other "metallic" DMS materials 
where the ferromagnetic interaction between the impu- 
rity local moments is mediated by itinerant carriers (elec- 
trons or holes). Our theory is a suitable 2D generalization 
of both the continuum lattice version of the highly suc- 
cessful RKKY mean field theory (MFT) that has earlier 
been applied [1] to 3D metallic Gai_ x Mn x As systems. 
One of our important findings is that the ferromagnetic 



transition ("Curie") temperature for the 2D DMS sys- 
tems typically tends to be substantially less than the cor- 
responding 3D case with equivalent system parameters. 
In particular, T c in the 2D case is found to be compara- 
ble (T c ~ Tp) to the Fermi temperature (Tp = Ep/ks, 
where Ep is the 2D Fermi energy) of the 2D hole sys- 
tems. This necessitates that the full finite temperature 
form for the 2D RKKY interaction be used in calculat- 
ing DMS magnetic properties, further suppressing T c in 
the system. In fact, this general lowering of the 2D DMS 
T c compared with the corresponding 3D case is our cen- 
tral new theoretical result. This implies that spintronic 
applications involving 2D DMS heterostructures will be 
problematic since the typical T c (at least for the currently 
existing DMS materials) is likely to be far below room 
temperature (T c < 100K). A related equally important 
theoretical finding is that, although the continuum vir- 
tual crystal mean field approximation much used for 3D 
Gai^^Mn^As physics [4] would predict that the T c in 
2D DMS systems (being proportional to the 2D density 
of states) is independent of 2D carrier density, there is 
an intrinsic carrier density dependence of T c even in the 
strict 2D system arising from the density and temper- 
ature (i.e. T/Tp) dependence of the finite-temperature 
effective RKKY interaction. 

We specifically consider the so-called delta-doping 
schemes for localizing Mn dopants, the species respon- 
sible for DMS ferromagnetism, to very narrow layers [3] . 
It is this situation which we examine in this letter, spin- 
5/2 Mn impurities confined to a single plane in GaAs; 
although the experimental reality may well be one in 
which a significant number of impurities lie off this plane 
and/or the Mn occupy a different crystallographic plane, 
our results are qualitatively correct in these cases as well. 

Although the Hohcnberg-Mermin- Wagner theorem 
precludes long range order in 2D systems with Heisenberg 
spins, the theorem applies only for the case in which the 
coupling between spins is absolutely isotropic. As has 
been shown both formally and numerically [5-7], even 
a small amount of anisotropy is sufficient to stabilize 
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long range order at finite temperatures. We examine the 
impact of anisotropy explicitly by studying the classical 
anisotropic 2D Heisenberg model TL = —JoJ2{ij) [SfSj + 
SfSj + (7 + l)SfSj], where J is the ferromagnetic ex- 
change constant, 7 is the anisotropy parameter, and the 
spins occupy a 2D square lattice; the sum is over nearest- 
neighbors only. Using a variant of the Wolff cluster 
Monte Carlo technique (with modifications appropriate 
for a finite 7 [8]), we calculate Curie temperatures by 
finding the intersection of Binder cumulants (given by 
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I(|M| 4 )/(|M| 2 ) 



for two different system 



sizes (L — 40 and L = 80). The results are shown in 
Fig. 1, where it is readily evident from the main graph 
and the inset (depicting T c for a smaller range of 7 val- 
ues) that even a small amount of anisotropy yields Curie 
temperatures in the vicinity of Jo/fcs, a value in reason- 
able agreement with T^f FT = | Jo/fcs, the corresponding 
mean field result for the nearest neighbor 2D Heisenberg 
model. Previous theoretical studies also examined ferro- 
magnetism in the 2D DMS context [10,11], without, how- 
ever, considering the Hohenberg-Mermin- Wagner theo- 
rem or the anisotropy issue. 




FIG. 1. Curie Temperatures plotted versus 7, the 
anisotropy parameter. The inset displays T c values for a much 
smaller range of anisotropies. In both images, the Monte 
Carlo errors are smaller than the graph symbols. The black 
line is included only as a guide to the eye in both the main 
graph as well as the inset. 

We apply the lattice MFT developed previously [12], 
but we also find intriguing results in the context of sim- 
ple continuum Weiss MFT which demonstrate that it is 
essential to incorporate finite temperature effects in the 
effective carrier-mediated interaction between Mn impu- 
rity moments. We assume the 2D hole gas to be confined 
in the same plane as the Mn dopants- it is straightfor- 
ward to consider [10,11] a spatial separation between the 
dopants and the holes as well as to consider the quasi-2D 
confinement for the holes. These additional complica- 



tions would lower T c below the strict 2D limit considered 
in our model. 

Our theory is constructed for two dimensional DMS 
systems in the metallic limit with itinerant carriers (we 
assume the carrier-mediated effective Mn-Mn indirect ex- 
change interaction to be of the RKKY form). However, 
including a finite carrier mean free path in our theory 
allows us to take into account the dependence of the 
magnetic behavior of our system on the carrier transport 
properties. In fact, using an exponential cutoff in the 
range of the RKKY function permits us to treat the long 
and short-range magnetic interaction regimes simply by 
varying the cutoff parameter I. In our system, salient pa- 
rameters include the Mn local moment concentration (x), 
the free carrier density (n c ), and the exponential cutoff 
scale I associated with the carrier mean free path. The 
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2D Fermi temperature Tp — ~ — depends on n c since 
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the 2D kp oc n c , where kp is the Fermi wave vector. 
We take the carrier effective mass to be m* = 0.4m e for 
all our numerical results. All quantities (Tp,kp,r, etc.) 
appearing in our theory refer to 2D parameters. 
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FIG. 2. Temperature dependent 2D RKKY range func- 
tion xi^Fr, T/Tf) as a function of Uft for several temper- 
atures T/Tf = 0.0,0.5, and 1.0 (corresponding to the solid, 
dashed, and dot dashed lines respectively). Inset (a) also 
shows x(kFr,T /Tf), and inset (b) portrays the effective fi- 
nite temperature coupling J g-. 

Our effective Hamiltonian describes the Mn-Mn mag- 
netic interaction between classical Heisenberg spins Si on 
a 2D lattice: 



(i) 



where Si is the ith Mn local moment of spin 5/2. In our 
lattice MFT we assume the Mn dopants to lie entirely on 
the [100] plane of the GaAs zinc-blende crystal lattice, 
with the impurities occupying the Ga sites, which form 
a square lattice with lattice constant a. 
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In our Hamiltonian, the carrier mediated RKKY in- 
direct exchange interaction describes the effective mag- 
netic interaction between Mn local moments induced 



by the free carrier spin polarization. J, 



RKKY 



(T) 



Jox(^-F r ! T/Tp), where x(kFf, T/Tp) is the temperature 
dependent range function which is obtained from the 



Tak- 



2D spin susceptibility, and Jo = 

ing J p< i — 0.15 eVnm 3 [13] and square lattice constant 
a = .4nm, we find Jo = 0.10 eV; we use this value for Jo 
throughout this work. 

For T = the 2D RKKY interaction is known exactly 
[14]: 

J* KKY = Jo(kpa) 2 [Jo(kFr)No{k F r) + J^kp^N^kpr)] 

(2) 

where J n (x) are the Bessel functions of the first kind, 
and N n {x) are Bessel functions of the second kind. The 
asymptotic form of the RKKY interaction for kpr ^> 1 
is J HKKY (r) oc Jo sa\{2kpr) / r 2 . At finite tempera- 
tures (i.e., for T/Tp ^ 0), it is not possible to obtain 
the RKKY range function analytically. Thus, we cal- 
culate the finite temperature RKKY interaction numer- 
ically. In Fig. 2, we show the temperature dependent 
range function x(kpr,T/Tp) as a function of kpr for 
various T/Tp values. One sees increasingly severe ther- 
mal damping of the RKKY oscillations with increasing 
T/Tp. Inset (a) of Fig. 2 provides a closer view of the 
RKKY oscillations depicted in the main graph, while in- 
set (b) displays the effective coupling constant, given by 
"^eff = / Jij KKY i r )dr. One finds, as expected, that J Q g 
decreases with increasing temperature due to the damp- 
ing of the range function. Since the magnetic properties 
of the 2D system are directly dependent on J e g, a finite 
T/Tp can have an important impact. 

In the lattice MFT (where we have in mind the previ- 
ously mentioned square lattice with lattice constant a), 
we calculate the Curie temperature with [12] 
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xJ2NiJ(n,T c ), 



where Ni and are the numbers and distances of the 
ith nearest neighbors, respectively. The continuum limit, 
which we examine first, is attained for I, kp 1 a, where 
a is the lattice spacing. We examine the large I limit 
(/ S> fc^ 1 ) and find that Eq. 3 becomes 



= 35-kxJq 
c ~ 6k B 
35irxJo 
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kp / x(k F rX/T F )rdr 
Jo 



9(t:/t f ), 
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where T* is the Curie temperature in the continuum 
limit, and g(T* /Tp) depends only on the ratio of T* to 



the Fermi temperature Tp. For T* < Tp, g(T*/T F ) -> 
1, and the dependence on carrier concentration is lost, 
and T* depends only on the impurity concentration 
x. (For convenience, we identify this temperature as 
T c °° = T c (Tp -» oo) = This result, with no de- 

pendence on n c , is proportional to the density of states 
at the Fermi Level Ep in two dimensions, and reflects 
the peculiarity of the 2D density of states being indepen- 
dent of carrier density. (The analogous continuum 3D 
result, also proportional to the density of states at Ef, 

1 /3 

is T* oc xn c , and does depend on the carrier density). 
For the more general case of finite T* /Tf, one can read- 
ily calculate the Curie temperature self-consistcntly via 
numerical means. Fig. 3 shows T* curves as a function 
of carrier density for various values of x, the Mn con- 
centration. A salient feature of the curves is a marked 
dependence on n c , in contrast to what one finds when the 
(T=0) RKKY formula is used (where there is no depen- 
dence at all on n c ). Eventually, the curves saturate as n c 
is increased. This is natural, since Tp increases with n c , 
and ultimately T* /Tp approaches zero. This Tp — > oo 
limit of the continuum MFT is the actual virtual approx- 
imation limit. 
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FIG. 3. 2D Curie temperature curves in the continuum 
MFT for various values of the magnetic impurity density nt 
(m = xaT 2 ). Here is the continuum T c calculated for 



(3) Tf — oo, i.e. the virtual crystal approximation result. 



Returning to the lattice MFT, which explicitly takes 
into account the discreteness of the square lattice, we 
also incorporate the effects of the finite carrier mean free 
path I by including an exponential cutoff in the range of 
the RKKY interaction. T c curves are shown in Fig. 4 
for several Mn dopant concentrations. As in Fig. 3, the 
Curie temperature increases monotonically in n c over 
the experimentally accessible range of carrier densities. 
However, for considerably higher n c (i.e. approaching 
10 14 cm -2 ), nonmonotonic behavior is seen in T c . This 
is evident in the inset of Fig. 4, where the Curie temper- 
ature curves seem to achieve saturation for intermediate 
carrier densities and ultimately begins to decrease as the 
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length scale kp 1 of the RKKY oscillations shrinks relative 
to the lattice constant a. Eventually, for kp 1 ~ a, the dis- 
crete nature of the lattice sum in Eq. 3 has a strong effect, 
leading to a considerably smaller T c than the continuum 
value, T*. This result for 2D T c including disorder and 
finite temperature RKKY effects is one of our main new 
results. 



Curie Temperature in Eq. 3 with the continuum version, 
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FIG. 4. Lattice MFT 2D Curie temperature curves for var- 
ious values of :r. In the inset, T c curves are shown for a much 
greater range of carrier density n c . In both the main graph 
and the inset, l/a = 5. 

T c is strongly affected by the finite size of I; in fact, 
one sees that for each of the carrier concentrations rep- 
resented in Fig. 5, T c initially increases sharply with 
increasing I, eventually saturating for i/n > 1. It is 
informative to examine the regime a -C l,kp , where 
one can operate in the continuum limit. We introduce 
T] = kpl as an important dimensionless variable; as will 
be seen, the ratio T* /Tp tends to zero as rj becomes 
small. For T c <C Tp and kpr <C 1, a reasonable approx- 
imation for the finite temperature RKKY range func- 
tion is J RKKY {r, T) « J RKKY (r, 0)[(1 - a 2 r 2 - a 3 T 3 ) + 
(b 2 r 2 + b 3 T 3 ){k F r)}e- r / 1 where r = T /Tp, a 2 = 0.2153, 
a 3 = -0.140, b 2 = -1.333, and b 3 = 0.862 (The a, and 
bi have been calculated numerically). Since J RKKY (r, 0) 
varies slowly with kpr for small kp, we can replace the 
Bessel functions in Eq. 2 with a series expansions in kpr, 
and the result for small values of this expansion variable 
is 

J RKKY (r,T*) « J k 2 p/n [ai + a 2 (k F r) 2 ] 

x [1 - (a 2 r 2 + a 3 r 3 ) + (b 2 r 2 + b 3 T 3 )(kpr)] e~ r l\ (5) 

where ai ee [-1/2 + 7 + ln(fc F r/2)], and a 2 = [-3/16 + 
7/4+1/4 ln(fc F r/2)]; 7 = .57722 ... is the Euler constant. 
In calculating T* (in the continuum case), we assume 
that although I <g; kp 1 , I » a. This additional condition 
allows the replacement of the discrete formula for the 



and one has T* 
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out the integration and solving for T*, we find 



T* = Jo 
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where 77 = kpl <C 1 and only terms up to fourth order 
in 77 are shown. The expression for T* given in Eq. 6 
reveals a strong suppression of T* due to localization 
effects; when kp is small (i.e. in the small n c limit) the 
RKKY range function is very extended relative to a and 
I. As a consequence, most of the RKKY interaction is 
truncated by the exponential cutoff associated with the 
finite mean free path I. This truncation effect is so severe 
that T* is very small in comparison with Tp; to fourth 
order in r\ there are no corrections to arising from the 
finiteness of T* /Tp. Note that in the opposite limit of 
kpl ^> 1, which rarely applies to DMS systems which are 
at best "bad metals" , one can obtain the simple formula 
T c (k F l > 1) w T*[l - (kpl)^ 1 /n] by simply considering 
the disorder induced suppression of the 2D density of 
states. 
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FIG. 5. Lattice MFT Curie temperature curves for various 
values of n c plotted as a function of the mean free path l/a. 
In each case the same Mn concentration, x = .05, is used. 

In addition to T c , one can also use lattice MFT to cal- 
culate the magnetization M(T) [12]. The magnetization 
behavior is influenced by the Mn impurity concentration 
as well as the form of the effective interaction between 
Mn local moments, in principle specified by the carrier 
mean free path I and n c . For convenience, we study 
M(T/T C ); using the normalized temperature scale T/T c 
allows a systematic comparison of magnetization profiles 
corresponding to different values of the parameters I, n c , 
and x. An important trait of the magnetization profile is 
its degree of concavity; concavity in M(T) is a hallmark 
of an insulating system, while a convex profiles occur 
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well within the metallic regime [4] . Linear magnetization 
curves correspond to intermediate impurity densities and 
mean free paths. In terms of the magnetization, the con- 
cavity 7 is given by 7 = J t ^ M (T)dT, or the difference 
in the slopes of M(T). The sign of 7 indicates whether 
M(T) is convex (negative 7), concave (positive 7), or lin- 
ear (if 7 « 0). The temperatures t\ and ti are selected 
to encompass an intermediate temperature range, neither 
too close to T c nor to zero. A significant feature of the 
concavity plots shown in Fig. 6 is weak dependence of 
the concavity of the magnetization profiles on the carrier 
density; though the three values of n c range over two or- 
ders of magnitude (10 10 — 10 12 ) cm~ 2 , the 7 graphs lie 
very close to one another. One can also see that the pa- 
rameter range over which the M(T) profile is concave is 
much more extensive than predicted by lattice MFT for 
three dimensional DMS systems [12]. This is primarily a 
consequence of the geometry of the two dimensional lat- 
tice. The inset of Fig. 6 displays representative cases of 
the magnetization M(T). 
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FIG. 6. Concavity curves versus Mn concentration x for 
several carrier densities n c with l/a — 2.0. The inset shows 
representative examples of the magnetization, M(T). All 
M(T) curves are calculated for n c = 2.5 x 10 11 cm -2 . 



free path. In general, the 2D DMS T c is strongly sup- 
pressed compared with the corresponding 3D DMS T c , 
which is somewhat discouraging for spintronic applica- 
tions. 
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In conclusion, we have considered diluted magnetic 
semiconductors in two dimensions. We find that long- 
range ferromagnetic order can be stabilized by even a 
small amount of anisotropy invariably present in real- 
ity. We have found that, even at the level of continuum 
MFT, finite temperature effects in the carrier-mediated 
effective interaction between Mn impurity moments in- 
troduce a strong dependence on the density of carriers, 
where naive use of the zero temperature RKKY formula 
yields a result independent of n c . To take into account 
the discreteness of the strong positional disorder of the 
2D DMS system, we have employed a lattice MFT. Our 
lattice theory also provides a convenient framework for 
the inclusion of important physics such as the finite mean 
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